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FOREWORD 


NASA experience has indicated a need for uniform criteria for the design of space vehicles. 
Accordingly, criteria have been developed in the following areas of technology: 

Environment 

Structure 

Guidance and Control 
Chemical Propulsion 

Individual components of this work are issued as separate monographs as soon as they are 
completed. A list of the monographs pubUshed in this series can be found on the last page. 

These monographs are to be regarded as guides to design and not as NASA requirements, ex- 
cept as may be specified in formal project specifications. It is expected, however, that the 
monographs will be used to develop requirements for specific projects and be cited as the ap- 
plicable documents in mission studies, or in contracts for the design and development of 
space vehicle systems. 

This monograph replaces a monograph of May 1969 on the Earth’s upper atmosphere which 
provided a model for predicting atmospheric parameters at altitudes between 120 and 1000 
km. The current model given herein covers a wider range — 90 to 2500 km — and provides 
improved predictions of atmospheric parameters. 

This monograph was prepared under the cognizance of the Goddard Space Flight Center 
with Scott A. Mills of Goddard Space Flight Center and Robert E. Smith of Marshall Space 
Flight Center as program coordinators. The principal author was George Y. Lou of 
Northrop Services, Inc. The author is indebted to his colleagues, Richard L. King and 
Wallace W. Youngblood, for their valuable su^estions and criticism in the preparation of 
this monograph. 

Comments concerning the technical content of these monographs will be welcomed by the 
National Aeronautics and Space Administration, Goddard Space Flight Center, System Re- 
liability Directorate, Greenbelt, Maryland 20771. 
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MODELS OF EARTH’S ATMOSPHERE 
(90 to 2500 km) 

1. INTRODUCTION 


Atmospheric conditions encountered by a spacecraft in orbit about the Earth are important 
factors in space vehicle design, mission planning, and mission operations. Density is the pri- 
mary atmospheric property that affects the spacecraft’s orbital altitude, lifetime, and motion. 
Density directly affects the torques which result from aerodynamic interaction between the 
space vehicle and the atmosphere; such torques must be considered in design of spacecraft 
attitude control systems. Density scale height is required in heating calculations for space 
vehicles re-entering the Earth’s upper atmosphere. Density as well as chemical composition 
and temperature are needed in calculating a spacecraft’s drag coefficient. Chemical compo- 
sition and temperature also are required in the design of some experiment sensors to be flown 
in the upper atmosphere. 

Because of variability of atmospheric conditions with spatial location and solar activity, in- 
variant models of the Earth’s atmosphere (90 to 2500 km) would not be useful for most 
engineering applications. Therefore this monograph presents a computerized version of a 
model atmosphere developed by Jacchia in 1970 (ref. 1) that incorporates variation with 
solar activity and spatial location. The computerized model can provide at selected times 
and locations estimates of density, composition, temperature, molecular mass, pressure scale 
height, and density scale height for altitudes between 90 and 2500 km. 

This monograph replaces a monograph on the upper atmosphere which was a computerized 
version of Jacchia’s 1964 model (ref. 2). The current model has a range from 90 to 2500 
km; the model in the earlier monograph extended from 120 to 1000 km. The assumed 
boundary conditions at 1 20 km in the earlier monograph led to misrepresentation of atmos- 
pheric conditions below 200 km. Additional data from rocket probes and low-altitude sat- 
ellites made possible the development of the improved model recommended herein. 

In addition to the computerized technique, this monograph gives a quick-look prediction 
method that may be used to estimate the density for any time and spatial location without 
using a computer. Sample problems illustrate this method. 

Information in this monograph applies to altitudes between 90 and 2500 km; another mono- 
graph (NASA SP-8084) gives extreme atmospheric conditions at launch sites and another 
monograph is to give in-flight atmospheric extremes from the launch sites to an altitude of 
90 km. Other design criteria monographs on the Moon, other planets, and space vehicle 
technology are listed at the end of this monograph. 
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2. STATE OF THE ART 


Model atmospheres of the upper atmosphere have been developed that use assumed physical 
relationships in conjunction with observed phenomena (refs. 3 and 4). However, incomplete 
knowledge of the physics of the upper atmosphere limits the value of the resulting models 
for engineering purposes. 

This monograph concerns itself with model atmospheres that have empirical coefficients (de- 
' rived mainly from analysis of changes in periods of orbiting satellites) which relate neutral 
densities to solar flux, geomagnetic activity, time of day, latitude, and season. 

Since the publication of reference 5, a more detailed picture of the structure of the Earth’s 
upper atmosphere has emerged from additional satellite and rocket measurements, particu- 
larly in the altitude region below 200 km. As a result, the representativeness of the atmos- 
pheric models has been greatly improved. For example, the models recently published in 
references 1 and 6 have removed the long-recognized drawback of the invariance of the 
boundary conditions at 120 km; consequently, the atmospheric variability at that altitude 
and the altitude region immediately above it are represented better. In addition, significant 
improvements have been made by increasing the O/O 2 ratio at 1 20 km and representing the 
semiannual effect as a change in density rather than in exospheric temperature. From the 
space vehicle design and planning standpoint, the models of the Earth’s upper atmosphere 
given in references 1 and 6 appear to be the best available to date. Because the model given 
in reference 1 has been proven to yield better results in orbital lifetime predictions than its 
later version given in reference 6, the prediction method given in reference 1 was adopted in 
the present monograph. 

An empirical model recently has become available that gives more detailed information on 
densities of individual gas species on the basis of OGO-6 data (ref. 7). 


2.1 Earth’s Upper Atmosphere 

2.1.1 Terms of Reference 

For purposes of reference, the Earth’s atmosphere is often divided into regions. In the most 
commonly used system, these regions are differentiated by thermal stratification. The low- 
est region, extending from ground level to the first temperature minimum, is called the tro- 
posphere. The next region, which extends up to the level with maximum temperature, is 
called the stratosphere. The third region, which extends up to the second temperature mini- 
mum, is the mesosphere. The region directly above the mesosphere is the thermosphere. Be- 
yond the thermosphere is another region called the exosphere, but the latter is not defined 
by its temperature distribution. These terms and other terms that identify regions character- 
ized by other phenomena are given in the Glossary (appendix C). 


2 


Besides the foregoing divisions, the atmosphere often is divided into a lower and an upper 
region with the choice of altitude to separate them being somewhat arbitrary. In the present 
monograph, the term “upper atmosphere” is adopted to designate the portion of the atmos- 
phere above 90 km in altitude. 


2.1.2 Variation of Composition and 
Temperature with Altitude 

A brief description of the phenomena associated with the variations in composition and tem- 
perature with altitude is helpful for an understanding of the upper atmosphere. In the 
homosphere, the mixture of gasses is distributed vertically in quasi-hydrostatic equilibrium, 
with turbulent mixing keeping the relative composition of the atmosphere constant from the 
Earth’s surface to an altitude of about 90 km. Thus a given temperature profile in the 
homosphere can uniquely determine the corresponding density distribution to be associated 
with a pressure value at a given lower altitude. 

As the altitude increases, photodissociation processes cause the composition to change, re- 
sulting in a decrease of the molecular mass with height, primarily from replacement of mo- 
lecular oxygen by atomic oxygen. At higher altitudes, diffusive separation becomes domi- 
nant so that each individual gas is distributed by altitude according to its own molecular 
mass. Consequently, the abundance of lighter gasses, such as atomic oxygen, atomic nitro- 
gen, atomic hydrogen, and helium, decreases less rapidly with height than is the case for the 
heavier gasses. Above certain altitude levels these lighter gasses predominate over the heavier 
ones. It is then difficult to determine the distribution of mean molecular mass. Thus, a 
given temperature profile will not define a unique density profile in the heterosphere. 

Because of the long and short term variations in the upper atmosphere, there is no unique 
altitude that defines the base of the heterosphere. It occurs somewhere between 80 and 100 
km above the Earth and closely coincides with the base of the thermosphere. The thermos- 
phere extends to an average altitude of about 500 km. The upper limit of the thermosphere, 
i.e., the thermopause, varies with the level of solar activity. During periods of minimum 
solar activity, the thermopause occurs at about 350 km; during periods of maximum solar • 
activity, it can be as high as 700 km. Above the thermopause is the region known as the ex- 
osphere. In this region, atoms of hydrogen and helium can break loose from the Earth’s 
gravitational field because of their lightness and energy and escape into outer space. In the 
lower exosphere, the most abundant atom is oxygen. At higher altitudes, the most abundant 
atom is helium, and at still higher altitudes, atomic hydrogen. 

A basic feature of the structure of the thermosphere is the very steep increase in tempera- 
ture with height from 90 to 200 km. This is because of absorption of solar radiation 
primarily with wavelengths up to 2000 A, particularly in the extreme ultraviolet (EUV) 
portion of the spectrum (ref. 8). The absorbed radiation causes dissociation and ionization 
of particles, and consequent release of heat. For altitudes above approximately 200 km, the 
heat deposited in the atmosphere decreases with altitude and the atmosphere approaches 
isothermality if there is not appreciable heat input from magnetosphere. 
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2.2 Variations in Atmospheric Parameters 
With Location and Solar Condition 

The vast number of determinations of the drag exerted by the atmosphere on satellites and 
measurements made by density gauges and mass spectrometers aboard rockets and satellites 
have revealed seven different effects other than altitude that result in variations of density, 
temperature, and composition of the upper atmosphere. These observed variations are; 

( • Variations with solar activity 

• Diurnal variation 

• Variations with geomagnetic activity 

• Semiannual variation 

• Seasonal-latitudinal variations of the lower thermosphere 

• Seasonal-latitudinal variations of helium 

• Rapid density fluctuations probably associated with tidal and gravity waves 

Each of these variations have been taken into account in the model adopted herein except 
for the variations associated with tidal and gravity waves. More information can be found in 
references 1 , 6, and 9. 


2.2.1 Variation with the Solar Activity 

The solar ultraviolet radiation impinges upon the Earth’s upper atmosphere continuously, 
but with variable spectral intensity distribution according to solar conditions. It is the UV 
and EUV radiation that heats and causes structural changes in the upper atmosphere. There 
are two components of this radiation that generate heat energy. The one that relates to ac- 
tive regions on the solar disk varies from day-to-day, whereas the one that relates to the solar 
disk itself varies more slowly with the 1 1-year solar cycle. The atmosphere has been observed 
to react to each of these two components in a different manner (ref. 1). Separate values of 
the two components of the solar flux are not readily available; however, the observations 
made by OSO 1 in 1962 have confirmed that the integrated EUV flux and the 10.7-cm radi- 
ation (a quantity readily measurable at the Earth’s surface) are proportional to each other. 

In the thermosphere, the density is strongly influenced by the changing levels of solar activ- 
ity. The resulting density variations, in turn, affect the frictional air drag on satellites. The 
relationship between solar activity and the air drag exerted on satellites was first recognized 
by Jacchia in 1958 (ref. 10). From the study of the orbits of Sputnik III and Vanguard 1, he 
found that the air drag increased considerably in periods of high solar activity. During these 
periods of solar flares and other solar disturbances, the Sun’s ultraviolet and corpuscular radi- 
ation increase markedly; the result is greater heating of the Earth’s upper atmosphere. Sub- 
sequent studies of satellite drag data have verified a pronounced correlation between solar 
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activity and density variations (refs. 11 through 14). However, the reaction of the atmos- 
phere to variations of the Sun is not instantaneous. From analysis of the time lag between 
EUV flux from OSO-1 and exospheric temperatures computed by Jacchia and Slowey 
(ref. 15), Bourdeau et al. (ref. 16) concluded that the atmospheric response time was on the 
order of one day. Recently, Roemer (refs. 17 and 18) found a time lag of 1,0 ±0.12 days in 
the 355 to 710 km altitude range on the basis of data from six satellites. 


2.2.2 The Diurnal Variation 

1 

2.2.2.1 Density 

The diurnal variation of density in the Earth’s upper atmosphere has been observed by satel- 
lite drag measurements. From approximate equality between day and night densities at 200 
km, the excess density of day over night increases with altitude; at 600 km the daytime den- 
sity can be higher by a factor of eight (ref. 19). This day-to-night ratio has also been ob- 
served to be strongly latitudinal dependent; the larger values are found at lower latitudes and 
smaller values at higher latitudes (ref. 20). Results of satellite drag measurements have clearly 
shown that the density has a maximum around 1400 local solar time at a latitude approxi- 
mately equal to that of the subsolar point, and a minimum around 0300 hours at about the 
same latitude in the opposite hemisphere. The effect is believed to be caused by the absorp- 
tion of EUV radiation and by heat conduction of the neutral gas. The energy that is con- 
ducted downward into the lower thermosphere and mesosphere is mostly lost by radiation 
processes. At an altitude of about 1 20 km, the time constant for heat loss by conduction is 
on the order of one day or greater (ref. 1 9). Thus, the diurnal variation is not a predominant 
phenomenon at lower altitudes. 

Theoretically, a satisfactory explanation of the diurnal process can be obtained by solving 
the hydrodynamic and thermodynamic equations; however, attempts made by Harris and 
Priester (ref. 21) and later by Mahoney (ref. 22) gave results that are not in agreement with 
the observed phenomena; e.g., the calculated diurnal bulge occurs about 1700 local solar 
time compared with the observed occurrence at about 1 400. This discrepancy has been ex-, 
plained to be the result of not taking into account the horizontal energy transport. Subse- 
quent studies especially those by Volland, Stubbe, and Mayr (refs. 3 and 23 through 26), on 
the effect of horizontal wind on the shape of the diurnal temperature and density distribu- 
tion did indicate a reduction of the time lag between theoretical and observed bulges, but 
the problem has not been totally resolved. 

Also, related studies have been made on derivation of thermospheric winds from pressure 
gradients of the diurnal bulges (refs. 27 through 33). Harris (ref. 29) has shown that the 
effect of ion drag on neutral particles has an important bearing on wind velocities. Use by 
Harris (ref. 29) of the electron density profile of Chandra (ref. 34) and the collision fre- 
quency between neutral and ions of Chapman (ref. 35) led to a shift of the calculated bulge 
one hour earlier to 1 600 local solar time. 
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2.2.2. 2 Temperature 


The diurnal variation of temperature in the thermosphere has been obtained from analyses 
of observations of satellite drag and the ground-based incoherent scatter technique. There is 
a large discrepancy in the phase of the variations as inferred from these two types of obser- 
vations. Chandra and Stubbe (ref. 36) concluded from a case study that this discrepancy can 
be fully reconciled by introducing temperature and density perturbations at the base of the 
thermosphere; however, their analysis was later re-examined by Cummack and Butler (ref. 
37) who concluded that the reaction of the thermosphere to temperature variations at the 
base of the thermosphere is small. The discrepancy has not been reconciled. 

f 

The amplitude of the diurnal temperature variation in the thermosphere is not constant. 
Through analyses of satellite drag data, Jacchia and Slowey (ref. 38) have shown that the 
maximum-to-minimum ratio of exospheric temperatures was about 1.32 from 1958 to 
1963, then it dropped to 1.26 in the middle of 1963. This low value was maintained until 
1967. The ratio of day-to-night exospheric temperature was also found not to depend 
completely upon solar activity although it tends to be higher during maximum solar activity 
and lower during minimum solar activity. A recent study (ref. 1) has concluded that this 
ratio has a fairly good correlation with the yearly running mean of the geomagnetic activity 
index Kp , but is independent of latitude. The numerical value of the ratio varies between 
1.27 and 1.40 and has an average value of 1.31. 


2.2.3 Variations with Geomagnetic Activity 

2.2.3. 1 Origin 

Geomagnetic storms usually occur when clouds of charged particles collide with the Earth’s 
magnetosphere. These charged particles are believed to be ejected from the Sun during the 
course of a solar flare which is generally a short-lived phenomenon. As a result, a large 
amount of solar radiation is emitted by the flare region which subsequently heats the Earth’s 
atmosphere. The heating mechanism is not well-understood, but there have been some in- 
vestigations (refs. 39 through 41). 

2.2.3.2 Time Lag 

Satellite observations have shown a time lag between the onset of geomagnetic disturbances 
and the subseqeunt observed increases in density. Early studies of Explorer 9 satellite drag 
data (refs. 1 5 and 1 7) showed this time lag to be about 5 hours. A recent investigation of 
high-inclination satellite data (ref. 42), however, has indicated that the atmospheric pertur- 
bation lags about 5.8 hours at high latitudes and that the time lag tends to increase towards 
the equator to 7.2 hours. From an analysis of LOGACS* data (ref. 43), the atmos- 
pheric density is found to respond to geomagnetic storms almost immediately near the polar 
region and to respond more slowly with increasing distance from the polar region. An aver- 
age time lag is presently taken as 6.7 hours (ref. 1). There is also evidence that even the 

*Low-G Accelerometer Calibration System 
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smallest variation in the magnetic field of the Earth, such as those observed during magnet- 
ically quiet days, is also reflected in density and temperature variations in the upper atmos- 
phere (refs. 44 and 45). There is recent evidence from the OGO-6 mass spectrometer that 
the atmospheric response to magnetic activity is dominated by N 2 increases at high latitudes 
(ref. 46), and this information may be reflected in a subsequent revision of this monograph. 


2.2.4 The Semiannual Variation 

Semiannual variation of the atmospheric density was first discovered by Paetzold and 
Zschorner (refs. 47 and 48) from satellite drag measurements in the altitude range of 210 to' 
650 km. It was later confirmed by Cook and Scott (ref. 49) at an altitude of 1 130 km and 
by King-Hele and Kingston (ref. 50) at 190 km. The outstanding characteristics of the semi- 
annual density variation are a primary minimum in July, followed by a high maximum in 
October, and a secondary minimum in January, followed by a secondary maximum in April. 
More recent analyses of satellite drag data (ref. 51) show that the semiannual variation in 
density in the 150 to 180 km altitude region is similar to that at higher altitudes. It has also 
been found (ref. 52) that the semiannual effect varies considerably from solar cycle to solar 
cycle, both in magnitude and in altitude dependence. According to Jacchia (ref. 2), the am- 
plitude of this variation is quite large at sunspot maximum but decreases toward sunspot 
minimum. 

The cause and mechanism of the semiannual variation are not yet fully understood. Early 
hypotheses are considered in references 48 and 53. References 54, 55, and 56 are recent 
attempts to explain this phenomenon. 


2.2.5 Seasonal'Latitu dinal Variations 
of the Lower Thermosphere 

Recent low-altitude, high-inclination satellite and rocket probe measurements have shown 
that the lower thermosphere is subject to a large seasonal-latitudinal variation in temperature 
and a smaller variation in density. The amplitude of the density variation was found to in- 
crease very rapidly with height from 90 km up to a peak somewhere between 105 and 120 
km (refs. 57 and 58) and then to decrease with altitude to 200 km where no appreciable 
seasonal-latitudinal variation has been observed (refs. 38 and 59). At higher altitudes, how- 
ever, a density anomaly in the auroral region has been reported by Newton and Pelz 
(ref. 60), Philbrick and Mclsaac (ref. 61), and Allan (ref. 62). A density bulge at altitudes 
below 200 km in the auroral region had been previously found by Jacobs (ref. 63) and 
Marcos et al. (ref. 64); this has been confirmed by recent analyses of the OVI 15 satellite 
drag data (ref. 65). 


2.2.6 Seasonal-Latitudinal Variation of Helium 

Reber and Nicolet (ref. 66) first noted the enhancement of helium toward the winter 
hemisphere in the analysis of mass spectrometer data from Explorer 17. 


7 


From the analyses of orbital decay of the Explorer satellites, Keating and Prior (ref. 67) 
observed a winter density bulge in the polar exosphere which they attributed to the peak 
helium concentration at high latitudes in the winter hemisphere (refs. 68 and 69). This 
finding has been confirmed by other independent satellite drag measurements (refs. 38 and 
70) and by analyses of the intensity of the 10830 A emission line (refs. 71 through 73). 
Recent mass spectrometric studies (refs. 74 through 77) also have revealed that the concen- 
tration of helium in the lower thermosphere tends to increase in the winter at high latitudes. 

Through an analysis of the intensity of the 10830 A emission line, Bitterberg et al. (ref. 78) 
deduced that helium varied seasonally in one year by a factor of 3 to 4 above 500 km. They 
'pointed out that clear maxima in the helium concentration were observed repeatedly in the 
month of December and January and then were followed by a rather steep decline. Minima 
in the helium concentration were usually observed in April and May, but some early ones 
were observed in March. 

The formation of this helium bulge over the winter pole has been explained by a seasonal 
subsidence of the level at which the diffusion of helium begins (refs. 68 and 79). It has been 
shown (ref. 38) that a change of this level by 5 km could change the amount of helium in 
the thermosphere by a factor of 2. However, the mechanism of this winter helium bulge and 
its latitudinal dependence are still under investigation. 


2.2.7 Rapid Density Fluctuations Associated 
with Tidal and Gravity Waves 

Rapid density fluctuations propagating throughout the upper atmosphere have been detected 
by density gauges aboard the Explorer 32 satellite (ref. 80) and by accelerometers aboard the 
OVI-15 satellite (ref. 81) in the altitude range from 120 to 510 km. These fluctuations are 
believed to be caused by tidal and gravity waves. They appear to be more prevalent at higher 
latitudes near the auroral region in the early morning or late evening hours and are observed 
to propagate mainly in a North-South direction with maximum horizontal wavelengths on 
the order of 130 to 520 km (ref. 80). Their vertical wavelengths increase with hei^t. Am- 
bient density perturbations with a half-amplitude of 50 percent from smooth density values 
during a short duration have been observed (ref. 80). Theoretical studies of the density vari- 
ations associated with tidal and gravity waves have been conducted (refs. 82 through 89). 

Because of insufficient knowledge, tidal and gravity effects have not been taken into account 
in the atmospheric model recommended herein. 


2.3 Determination of Atmospheric 
Properties (90 to 2500 km) 

Semi-empirical models that use empirical coefficients to correlate neutral densities with solar 
flux, geomagnetic activity index, time of day, and other parameters are the best engineering 
tools at present for predicting the atmospheric conditions encountered by a spacecraft in 
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Earth orbit. Examples of this kind of model are given in references 1, 2, 5, 6, 48 and 90 
through 93, but those given in references 1 and 6 are the most recent ones and are the best 
models available now. Of these two, the model given in reference 1 has shown better results 
in orbital lifetime predictions than its revision given in reference 6. Therefore, the prediction 
method given in reference 1 is adopted in the present monograph. It should be noted, how- 
ever, that the model in reference 1 is derived mainly from satellite drag data which are aver- 
aged over a considerable range in time and space. Therefore, in comparisons of model pre- 
dictions with instantaneous measurements, some divergence is to be expected. 


2.3.1 Analytical Method 

The prediction method given in Jacchia’s 1970 model (ref. 1) basically is a static diffusion 
model which defines temperature and chemical composition and provides densities in agree- 
ment with satellite drag observations and, to a lesser degree, with rocket probe measurements 
from 90 to 1 100 km. In practice, densities are derived from the empirically determined tem- 
perature profile and assumed constant boundary conditions at 90 km. Mixing is assumed to 
prevail to an altitude of 105 km and any change in the mean molecular mass below this level 
is assumed to result only from dissociation of oxygen. 

The distribution of mean molecular mass between 90 and 105 km is determined empirically 
in such a way that it results in a ratio of atomic oxygen to molecular oxygen of 1.5 at 120 
km. The model also assumes that diffusive equilibrium takes place above 105 km. All of the 
recognized variations in the upper atmosphere as described in section 2.2 except the one as- 
sociated with tidal and gravity waves are included in the model. More information on the 
structure of the model is given in reference 1. 

A computational procedure for the prediction method is given in appendix A. The compu- 
tation has been computerized and the output of the computer program gives temperature, 
mass density, number density, molecular mass, pressure scale height, and density scale height 
as functions of altitude. 


2.3.2 Table Loolc-Up Method 

A table look-up density model is provided for cases when only an estimate of atmospheric 
density is needed. This method eliminates the use of a computer; only a few simple hand 
calculations are required for computing the exospheric temperature and possibly one correc- 
tion term for density if the altitude of interest is below 170 km. The accuracy of the result- 
ing density from this simplified method should be very close to the one obtained from the 
analytical method except at heights where the helium becomes the predominant constituent 
at altitudes above 600 km. 

The look-up method gives only the mass density as a function of exospheric temperature. 
Procedures for calculating the exospheric temperature and finding the density from this 
method are given in appendix B where sample problems are illustrated. 
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3. CRITERIA 


Prediction models of the Earth’s atmosphere from 90 to 2500 km for use in space vehicle 
design and mission planning should be obtained in accordance with section 3.1 or 3.2. 

Section 3.1 should be used for predictions of density and associated atmospheric parameters; 
section 3.2 should be used when a quick estimate of atmospheric density is needed. 

Both methods (sections 3.1 and 3.2) can be used to obtain either mean density models or 
models having reasonable upper density extremes. Mean density values are obtained if nom- 
inal values for predicted solar and geomagnetic activity are used for inputs; upper density 
values result from using plus 2o values for predicted solar and geomagnetic activity. Table 1 
gives sample predicted solar and geomagnetic data. 

The upper density model obtained from the plus 2a values, however, does not account for 
short term surges in geomagnetic activity, which are usually of 6 to 1 2 hours duration. When 
experiments or subsystems are considered to be sensitive to such short term effects, a geo- 
magnetic index of 400 should be used with the predicted plus 2a solar flux values to obtain 
an upper density model associated with extreme geomagnetic conditions. 


3.1 Method for Predicting Density and Associated 
Parameters of Earth’s Atmosphere (90 to 2500 1cm) 

The prediction method given in appendix A (based on Jacchia’s 1970 model) should be used 
to predict atmospheric density, temperature, chemical composition, molecular mass, pressure 
scale height, and density scale height for any time and spatial location. 

Symbols, inputs, and calculation procedures are given in appendix A. Copies of the com- 
puter program for this method and the required inputs of the predicted solar and geomag- 
netic activity (issued monthly) are available upon request from the NASA Marshall Space 
Flight Center, Code S&E-AERO-Y, Huntsville, Alabama 35812. 


3.2 Method for Estimating Density of Earth’s 
Atmosphere (90 to 2500 1cm) 

The Quick-Look Density Model given in appendix B can be used to obtain an estimate of at- 
mospheric density for any time and spatial location. Although the Quick-Look Density 
Model provides density only, other atmospheric parameters may be interpolated from the 
tables given in reference 1 by using the exospheric temperature that is calculated by the 
Quick-Look Density Model. 
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The computation procedure of the Quick-Look Density Model is based upon the physical re- 
lationships given in appendix A, but all of the equations have been replaced by tables to 
eliminate the need for a computer. If interpolations are made with a three decimal accuracy, 
the density obtained from the Quick-Look method will be within 2 percent of that obtained 
by the method given in appendix A except at high altitudes above approximately 600 km 
where helium is the predominant constituent of the neutral gas. In this case, the prediction 
method given in appendix A is recommended. Otherwise, the reader is referred to reference 
1 in which a simple method can be found to estimate density values at high altitudes with a 
correction term for helium variations. 

Symbols, inputs, calculation procedures, and sample problems are given in appendix B. 
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TABLE 1 


AN EXAMPLE OF NOMINAL AND 2 a PREDICTIONS OF SUNSPOT NUMBERS, 
MEAN 10.7 SOLAR FLUX, AND GEOMAGNETIC INDEX* 


(Calculation by MSFC Solax Prediction Program) 


Calendar Year 

Sunspot Number 

Mean 10.7 cm Solar Flux 


Nominal 

+ 2 a 

Nominal 

+2a 

Nominal 

+ 2 o 

19 G 8. 25 

107.30 

119.34 

153.76 

165.40 

13 

23 

19 C 8. 50 

107.49 

130. 37 

153.96 

176.07 

13 

23 

1968.75 

107.74 

140.40 

154. 18 

185. 77 

13 

23 

1969.00 

104.78 

143. 12 

151. 33 

188.40 

13 


1969.25 

103.31 

147.75 

149. 90 

192. 88 

13 


1969. 50 

102.57 

146. 09 

149.18 

191.27 

8 


1969.75 

77.64 

106. 36 

125.07 

152. 85 

8 

23 

1970.00 

94.67 

139.76 

141. 55 

185. 15 

8 

23 

1970.25 

90. 56 

128.30 

137.57 

174.06 

8 

23 

1970. 50 

82.42 

116. 91 

129.70 

163.05 

8 

23 

1970.75 

77.64 

106. 36 

125.07 

152.85 

8 

23 

1971. 00 

70.60 


118.27 

141.49 

8 

23 

1971.25 

62.48 


111. 55 

132. 54 

8 

23 

1971. 50 

57.16 


107. 15 

127.89 

8 

17 

1971.75 

51.17 

75.38 

102.21 

122. 39 

8 

17 

1972. 00 

46.04 


97. 98 

119.84 

8 

17 

1972.25 

42.73 


95.25 

118.43 

8 

17 

1972. 50 

38.23 


91. 54 

113.31 

8 

17 

1972.75 

34.02 

mam 


U 2.04 

8 

17 

1973.00 

31.43 

61.69 

85. 93 

110.00 

8 

17 

1973.25 

27.77 

57. 58 

84. 66 

107. 50 

8 

17 

1973. 50 

24.91 

53. 80 

82.95 

104. 39 

8 

17 

1973.75 

22.37 

50.18 

81.42 

M 01.40 

8 

17 

1974.00 

18.78 

44.82 

79.72 

96.98 

6 

17 

1974.25 

16.19 

41.14 

77.71 

93. 94 

6 

17 

1974. 50 

14.70 

33.74 

76. 82 

91.96 

6 

17 

1974.75 

12.31 

33.04 

75.39 

87.26 

6 

17 

1975.00 

10.94 

28. 22 

74. 56 

84. 93 

6 

17 

1975.25 

10.80 

24.65 

74.49 


6 

17 

1975. 50 

10.36 

21.66 

74. 22 


6 

17 

1975.75 

11.76 

26. 13 

75.06 


6 

17 

1976.00 

13.96 

33.64 

76. 38 

87.75 

6 

17 

1976.25 

17.38 

44.03 

78.43 

96. 33 

6 

17 

1976. 50 

23.24 

62.04 

81.95 

111.18 

8 

17 

1976.75 

30.81 

83. 81 

85, 42 

1 31 . 05 

8 

23 

1977.00 

39.26 

105.91 

92. 39 

152.41 

8 

23 

1977.25 

48.44 

123.83 

99.97 

169.74 

8 

23 

1977.50 

55.66 

135. 55 

105.92 

181.08 

8 

23 

1977. 75 

63.86 

149.23 

112. 68 

194.30 

8 

23 

1978.00 

73.32 

165.80 

120. 90 

210. 33 

8 

23 

1978.25 

81.35 

178. 53 

128.66 

222,64 

8 

23 

1978.50 

88.94 

187.01 

136.01 

230. 84 

13 

23 

1978.75 

90. 88 

186.67 

1.37.88 

230. 51 

13 

23 

1979.00 

91.12 

189.24 

138. 11 

232.99 

13 

23 


♦These data are inferred from predicted mean sunspot numbers by linear regression techniques (ref. 94) which were 
modified to make quarterly predictions (ref. 95) 
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APPENDIX A 

SUMMARY OF METHOD FOR 
COMPUTING MODEL ATMOSPHERES 

The computational procedure for obtaining predicted models of the Earth’s atmosphere for any time and 
spatial location is given below. A computer program that has been developed for this procedure is avail- 
able upon request from the NASA Marshall Space Flight Center Mail Code S&E-AERO-Y, Huntsville, 
Alabama 35812. 

Predictions of the 10.7-cm solar flux and geomagnetic activity are required as inputs. The nominal values 
of such data together with estimates of ± 2 a variations (95% confidence envelope) for each quarter, 10 
years into the future, are issued monthly and also should be obtained from the foregoing address. Table 1 
is an example of such data. The mean solar flux should be used twice, as an input for the daily (F) and 
81 -day mean (F) solar flux.* 


*In obtaining models of the Earth’s atmosphere for any time in the past, observed daily and calculated 81 -day mean solar 
flux values are used for F and F. 
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SYMBOLS 


A 

AV 

A2 


Parameter used in calculating diurnal correction for exospheric 
temperature (equation A-11) 

Avogadro's number, 6.02257 x 10^^ mole ^ 

Parameter used in calculating temperature for altitude greater than 
125 km (equation A-17) 

Geomagnetic activity index at 6.7 hours before computation time, 
unitless 


B. Coefficients for calculatlnv mean molecular mass ^eonarion A-1 81 . 

; _____ _ ^ ^ J __ ~ » r 

unitless 


DD Day number after Jan. 1 (input), days 

DDD Seasonal-latitudinal density correction in the lower thermosphere 
(equation A-20) , gm/cm3 

DENJ Mass density at geometric altitude Z with seasonal-latitudinal 
correction (equation A-21) , gm/cm3 

-9 

DENO Assumed mass density at lower boundary (90 km altitude), 3.46 x 10 
gm/ cm3 

DENS Mass density at geometric altitude Z without seasonal-latitudinal 
correction, gm/cm^ 

DS Declination of Sun (equation A-6), deg 
EM Atmospheric molecular mass, unitless 

-22 2 

F Daily 10.7-cm solar flux (input), 10 watts/m /cycles/second 

F 81-day mean of F, ending on date atmospheric property desired (input), 

10-22 watts/m3/cycles/sec 

FK Universal gas constant, 8.31432 Joule/mole-°K 

G Acceleration of gravity at geometric altitude Z (equation A- 19) , 

cm/sec2 

GP Greenwich meridian position (equation A-3) , deg 
HRA Hour angle of Sun (equation A-8) , deg 
i Denotes an individual constituent 

J Julian date (equation A-1) , days 

J* Computed parameter used in computing GP (equation A-2), years x 100 
k Boltzmann's constant, 1.380 x 10“^^ erg/°K 

LAX Latitude of computation point (input), deg 

LNG Longitude of computation point (input), deg 
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LS Celestial longitude (equation A-5), radians 
MM Greenwich mean time from 0000 GMT (input), minutes 

MO Mean molecular mass at 90 km altitude, unitless 

M(i) Molecular mass of constituent 1, unitless 

N(i) Number density of constituent i, cm 

N(TOT) Total number density (equation A- 31), cm"^ 

-3 

N'(He) Helium number density with seasonal-latitudinal correction, cm 

-3 

PAR Number of particles per unit volume, cm 

Q(i) Content of constituent i in fraction by volume, unitless 

R Parameter used in calculating diurnal correction for exospheric 

temperature (equation A-11) , unitless 

RAP Right ascension of computation point (equation A-A) , deg 

RAS Right ascension of Sun (equation A-7), deg 

g 

RE Radius of Earth, 6.356766 x 10 cm 

SH Atmospheric pressure scale height (equation A-34) , km 

SP Atmospheric density scale height (equation A-35) , km 

TAU Angle between computation point and density bulge (equation A-9), deg 

TC Nighttime minimum global exospheric temperature (equation A-10) , 

TE Exospheric temperature (equation A-14) , 

TG Geomagnetic activity correction for exospheric temperature (equation 
A-12), 

TL Diurnal correction for e^tospheric temperature (equation A-11), °K 

TS Semiannual correction for exospheric temperature (equation A-13) , ®K 

XX Temperature at inflection point, Z“125 km (equation A-15) , °K 

TZ Temperature at geometric altitude Z (equations A-16 and A-17) , ®K 

TO Assumed temperature at lower boundary 90 km altitude, 183 ®K 
T1 Parameter used in calculating temperature (equation A-16) , ®K/km 

3 

T3 Parameter used in calculating temperature (equation A-16) , ®K/km 

4 

T4 Parameter used in calculating temperature (equation A-16) , ®K/km 
W(i) Mass of constituent i (input), gm/mole 
Y Length of tropical year, 365.2422 days 

YR Year (input) , years 

n Parameter used in calculating the diurnal correction for exospheric 

temperature (equation A-11) , deg 

e Parameter used in calculating the diurnal correction for exospheric 

temperature (equation A-11), deg 
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T Parameter used in calculating semiannual correction for exospheric 

temperature (equation A-13) , unitless 

a Thermal diffusion coefficient for atmospheric constituent, oi=-“0.38 

for helium, a=0 for the others, unitless 

AHe Semiannual correction for helium number density (equation A-29), 
unit less 


INPUT 

a^ Geomagnetic activity index at 6.7 hours before computation time, unitless 
DATE Date, month/day/year 
DD Day number since Jan. 1, days 
F Daily 10.7 cm solar flux, used as unitless 

F 81 -day mean of F, ending on date atmospheric properties desired, 

used as unit less 

LAT Latitude of computation point, north (+) , south (-), deg 
LNG Longitude of computation point, east (+) , west (-), deg 

MM Greenwich Mean Time from 0000 GMT, minutes 

YR Year 

z Geometric altitude, km 
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COMPUTATIONAL PROCEDURE 

1. SUN'S DECLINATION AND HOUR ANGLE 

A. Julian date (days) 

J » 2441683 + (YR-1973) x 365 + DD 

where 

YR = year (input) 

DD = day number after January 1 (input) 

B. J* parameter (years x 100) 

J - 2415020 


(A-1) 


(A-2) 


C. Greenwich meridian position (deg) 

GP = 99.6909833 + 36000.76854 (J*) + 0.00038708 (J*)^ 

+ 0.25068447 (MM) (A-3) 

where 

MM = Greenwich mean time in minutes (input) 

GP must be between 0 and 360 deg. 

D. Right ascension of computation point (deg) 

RAP = GP + LNG (A-4) 

where 

LNG “ longitude of computation point (input) 

RAP must be between 0 and 360 degrees 

E. Celestial longitude (radians) 

LS = 0.017203 (J-2435839) +0.0335 sin [0.017203 (J-2435839)] 

- 1-^10 (A-5) 
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F. Declination of Sun (deg) 

DS = arc sin [sin (LS) sin (23.45*)] 

G. Right ascension of Sun (deg) 

RAS = arc sin 

H. LS must be converted to degrees 

I. Put RAS in quadrant of LS 

J. Compute hour angle (deg) 

HRA = RAP - RAS 


tan (DS) 
tan(23.45*) 


(A-6) 


(A-7) 


(A-8) 


II. TEMPERATURE COMPUTATION 
A. Exospheric temperature 

1. Angle between bulge and computation point (deg) 

TAU = HRA - 37* + 6* sin (HRA + 43*) (A-9) 

where 

TAU must be placed between +180 and -180 degrees 

2. Nighttime minimum global exospheric temperature (*K) 

TC = 383 + 3.32 F + 1.8(F-F) (A- 10) 

where 

—22 2 

F = Dally 10.7-cm solar flux (input), 10 watts/m /cycle/second 

F » 81-day mean of F, ending on date atmospheric property 
desired (input), 10”^^ watts/m^/cycle/second 

3. Diurnal correction (*K) 

TL = TC(1 + R sin^'^ 9) ^1 + A cos^'° (A-11) 
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where 


A = R 


2.5 ^ 2.5 ^ 

cos n - sin 9 

1 + R sin^*^ e 


n - j|lat-ds| 

e » -I-Ilat+dsI 

LAT = latitude of computation point (Input), deg 


DS = Sun’s declination (equation A-6), deg 
R » -0.19 + 0.25 logj^Q F(t - 400*^) 

F(t - 400*^) Is the value of F 400 days before the computation 
date 

4. Geomagnetic activity correction (*K) 

TG - 1.0 ap + 100 [1 - exp(-0.08 a^)] (A-12) 

where a^ Is the value 6.7 hours before the computation time 

5. Semiannual correction (®K) 

TS = 2.41 + F[0.349 + 0.206 8ln(360“T + 226.5*)] 

sln(720*x + 247.6*) (A-13) 


where 


1 + sin 


[S60* (1^) + 342.3*]1 


T - ^ + 0.1145 


Y « length of tropical year, 365.2422 days 

6. Compute exospheric temperature (*K) 

TE » XL + TG + TS 


2.16 




(A-14) 
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B. Temperature at inflection point, Z = 125 km, (*K) 
TX = 444.3807 + 0.02385 TE - 392.8292 


exp (-0.0021357 TE) 


C. Temperature at geometric altitude levels (°K) 

1. For altitude greater than 90 km but less than 125 km 
TZ = TX + TKZ-125) + T3(Z-125)^ + T4(Z-125)^ 


where 


T1 


- (^) 


T4 = 3 


T3 = - 


ftx - 183 - 2(T1) (1^)] 


T1 
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+ 4(T4) 


(B 


3(35)‘ 

2. For altitude greater than 125 km 

TZ = TX + (A2) tan"^ {T1(Z-125)[1 + 4.5 x 10“^ 


(Z-125)^'^]}] 


where 


A2 = 2 (TE-TX) 


III. NUMBER DENSITY, MASS DENSITY, AND MOLECULAR MASS COMPUTATIONS 
A. For altitude equal to or less than 105 km 

1. Mean molecular mass at altitude Z (unitless) 

EM = B^ + B2(Z-100) + B^CZ-lOO)^ + B^(Z-IOO)^ 

+ B.(Z-IOO)^ + B,(Z-100)^ + B,(Z-100)^ 

J D / 


(A-15) 


(A-16) 


(A-17) 


(A-18) 
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where 


= 28.15204 

B 2 = -0.085586 

= 1.2840 X 10 "^ 

B, = -1.0056 X 10 ”^ 

4 

B^ = -1.0210 X 10 "^ 

B, - 1.5044 X lO"^ 
o 

B^ = 9.9826 X lO"® 

2. Mass density at altitude Z before seasonal-latitudinal correction 
(gm/cm3) 


DENS-DENO (2|) (®) exp 

where 


"' (EM)G 

"(FK)TZ 


(A-19) 


DENO is the assumed density at 90 km altitude 
3.46 X 10“9 gm/cm^ 

TO Is the assumed temperature at 90 km altitude 183 ®K 

MO is the mean molecular mass at 90 km altitude 28.82678 

G is the gravity at altitude Z which can be calculated as 

_2 

G = 980.665 ^1 + (cm/sec^) 

3 

RE is the radius of the Earth, 6.356766 x 10 km 

FK is the universal gas constant, 8.31432 Joule/mole - ®K 

3 

3. Seasonal-latitudinal correction for density (gm/cm ) 


DDD = 0.02(Z-90) exp [-0.045(Z-90) ] 

sin^(LAT) sin (DD + 100) (A-20) 
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4. 


Mass density at altitude Z with seasonal-latitudinal correction 
(gm/cm3) 

DENJ = DENS X 10^°^ (A-21) 

-3 

5. Total number of particles per unit volume (cm ) 


PAR 


AV (DENJ) 
EM 


(A-22) 


where 

23 

AV is Avogadro’s number, 6*02257 x 10 mole 

-3 

6. Number density for molecular nitrogen, argon, and helium (cm ) 


■ “Xi) (i^) 


(A-23) 


where 


i denotes N 2 , Ar or He 

Q(i) is the content of N^, Ar or He (fraction by volume) 


Q(N2) = 0.78110 
Q(Ar) = 0.00934 
Q(He) = 0.00001289 

-3 

7. Molecular oxygen number density (cm ) 


N(02) = PAR 


EM 


28.96 


[1 + 0(02)]-:^ 


(A-24) 


where 


0 ( 02 ) * 0.20955 


-3 

8. Atomic oxygen number density (cm ) 


N(0) = 2 (PAR) 


1 - 


EM 

28.96 


(A-25) 


B. 


For altitude greater than 105 km 

_3 

1. Hydrogen number density at 500 km altitude (cm ) 
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where 


N(H) 5 oo - anti log {73.13 - 39.4 log (TE) + 5.5 [log (TE) ]^}(A-26) 


log denotes common logarithm 


2 . 


-3 

Hydrogen number density for altitude greater than 500 km (cm ) 
B(H) - (ff) e*p (- (A-27) 


where 


M(H) is the molecular mass of hydrogen, 1.00797 
k is Boltzmann's constant, 1.380 x 10 erg/®K 


3. Number density for molecular nitrogen, molecular and atomic 
oxygen, and helium (cm~3) 


N(l) - N(1),„5 .xp (- 


(A-28) 


where 


1 denotes N 2 , O 2 , 0 or He 

N(i)io 5 is the number density of each constituent at 105 km 
altitude 


4. 


a is the thermal diffusion coefficient, a=-0.38 for helium, 
and a«=0 for the other constituents 

Seasonal-latitudinal correction for helium number density (unitless) 


AHe = 0,5 + 1.8 


/ 23.45°-DSN 
V 47.5“ ) 





45“+DSj 


2.5 


47.5“ 


sin 


4 /jff lain 

U ~ 2 ; 


(A-29) 


-3. 


6 . 


Helium number density with seasonal-latitudinal correction (cm ) 

N'(He) = N(He) (AHe) (A- 30) 

-3 

Total number density (cm ) 

N(TOT) - N(H) + N(He) + N(N2) + N(02) + N(0) (A-31) 
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3 

7. Mass density (gm/cm ) 

DENJ = N(H) W(H) + N(He) W(HE) + N(N2) W(N2) 
+ N(02) W(02) + N(0) W(0) 

where 


W(i) denotes mass of constituent 1 In gm/mole 


W(H) = 1.6731 X 10 


-24 


-24 


W(He) = 6.6435 x 10 

W(N2) = 4.6496 x 10~^^ 

W(02) = 5.3104 X 10~^^ 

W(0) = 2.6552 X 10~^^ 

8. Molecular mass (unitless) 


_ (DENJ) (AV) 
^ “ N(TOT) 


(A- 32) 


(A-33) 


IV. PRESSURE SCALE HEIGHT (km) 


SH = 


8.31432 X 10~^ (TZ) 


(EM) (G) 

V. DENSITY SCALE HEIGHT (km) 

SP = 


8.31432 X 10-3 (tz)/(EM) (G) 
TZ 1 


, . 8.31432 d 

G dZ 

OUTPUT 


EM 


Temperature, “K 


Number density of N 2 , 0^, 0, He, and H, cm 

-3 

Total number density, cm 

3 

Mass density, gm/cm 


(A- 34) 


(A-35) 


Molecular mass , unitless 
Pressure scale height, km 
Density scale height, km 
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APPENDIX B 

OUlCK-LOOK DENSITY MODEL 
WITH SAMPLE PROBLEMS 


The computational procedure and sample problems for obtaining an estimate of 
atmospheric density for a given time and location are given below. Required 
inputs include predictions of the 10.7-cm solar flux and geomagnetic activity 
for the future time desired. Such data are issued monthly and current 
predictions should be obtained from the NASA Marshall Space Flight Center, 
Code S&E-AERO-Y, Huntsville, Alabama 35812. 


SYMBOLS 

a 

P 

DS 

F 

F 

F(t-400^) 

LAT 

P 

R 

S 

TC 

TE 

TG 

TL 

TS 


Geomagnetic activity index at 6*7 hours before computation time, 
unitless 

Declination of Sun, deg 

-22 2 

Daily 10.7-cm solar flux, 10 watts /m /cycles/second 

81-day mean of F, ending on date atmospheric property desired, 

10“22 watts/m^/cycle/second 

- -22 2 
Value of F 400 days before computation time, 10 watts/m / 

cycle/second 

Latitude of computation point, deg 

Parameter used for computing seasonal-latitudinal density 
correction in the lower thermosphere (Table B-4), unitless 

Parameter used in calculating diurnal correction for exospheric 
temperature (equation A-11), unitless 

Parameter used for computing seasonal-latitudinal density 
correction in the lower thermosphere (Table B-4), gm/cm^ 

Nighttime minimum global exospheric temperature (equation A-10), 
°K 

Exospheric temperature (equation A-14), 

Geomagnetic activity correction for exospheric temperature 
(equation A-12), ®K 

Diurnal correction for exospheric temperature (equation A-11) , 
Semiannual correction for exospheric temperature (equation A-13), 
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Geometric height, km 

3 

Mass density, gm/cm 

Parameter used for computing semiannual correction of exospheric 
temperature (Table B-3) , ®K 

Date, month/day/year (January 1, 1975) 

Daily 10.7-cm solar flux (74.56) 

ftl— rloTT tnck Q -n erila'T -Pl iiv (l lx 

y, , m ^ , 

Value of F 400 days before computation time (80.8) 

Geomagnetic activity index (6) 

Latitude of computation point (45®N) 

Local time (1000) 

Altitude of interest (130 and 400 km) 

COMPUTATIONAL PROCEDURE AND SAMPLE PROBLEMS 

The following examples illustrate the step-by-step procedure for obtaining 
an estimate of the atmospheric density over a point with geographic latitude 
45®N on January 1, 1975, at 1000 hours local time, for two heights: Z * 130 

km and Z = 400 km. From the sample data in table 1, the predicted nominal 
geomagnetic and mean solar flux indices are found to be 6 and 74.56, 
respectively. The geomagnetic index from table 1 is used as the value occurring 
6.7 hours before the computation time; and the mean solar flux is used twice, 
as the 81-day mean of 10.7-cm solar flux and the actual flux on the day 
before.** It should be noted, however, that table 1 is only a sample; the 
most recent predictions of the 10.7-cm solar flux and geomagnetic activity 
should be obtained for predicting atmospheric density. 


Z 

P 

6T 

s 

INPUT* 

DATE 

F 

p 

F(t-400'^) 

a 

P 

LAT 

LT 

Z 


♦Inputs for sample problems in parenthesis. 

**In obtaining models of the EarUi’s atmosphere for any time in the past, observed daily and calculated 81 -day mean solar 
flux values are used for F and F, and the observed value of geomagnetic activity index at 6.7 hours before the computa- 
tion time is used for 3p. 
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I. EXOSPHERIC TEMPERATURE COMPUTATION 

A. Nighttime minimum global exospheric temperature (equation A-10) 
TC = 383 + 3.32 F + 1.8(F-F) 

= 383 + 3.32 X 74.56 + 1.8(74.56-74.56) 

= 630.5 “K 


B. Exospheric temperature with diurnal correction 

1. From Figure B-1, the declination of the Sim on January 1 is 

DS = -23.45“ 

2. For LAT = +45® and local time « 1000, Table B-1 gives TL/TC * 1.116 
However, Table B-1 is computed for R*0.31. The actual R, 
according to equation A-11, is 

R » -0.19 + 0.25 log^Q F(t-400‘^) 

- -0.19 + 0.25 logj^Q (80.8) 

= 0.287 


where 


F(t-400 ) is obtained from Table 1. 

3. Therefore, the exospheric temperature with diurnal correction is 


TL - 630.5 X 1.116 

=• 630.5 X 1.033 
= 651.3 °K 


/0 .31 - 0. 
V 0.31 


287 


C. Geomagnetic activity correction 

From Table B-2, for a *6, we find 

P 

TG =■ +47 ®K 
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D. Semiannual correction 

From Table B-3, for January 1, we have 

6T = -11.6 °K 
s 

Therefore, the semiannual correction is 

= 0.7456 (-11.6) 

= - 8.6 

E. Exospheric temperature (equation A-4) 

TE = TL + TG + TS 
= 651.3 + 47 — 8.6 
= 689.7 °K 

II. ATMOSPHERIC DENSITY INTERPOLATIONS 

TE, which is the final exospheric temperature, should be used to enter 
Table B-5 to obtain the atmospheric density at any desired altitude. 

A. Density at Z=130 km 

1. From Table B-5, for TE ■= 689,7 “K, we have log^^^ p = -11.121 

2. For seasonal-latitudinal correction, we use Table B-4 where 

S = 0.132, P = 0.989, sin^(LAT) = 0.500 
Therefore, 

Alogj^QP = (S) (P) sin^ (LAT) 

= 0.065 

3. Final density at Z = 130 km 

log^gP = -11.121 + 0.065 
- -11.056 
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p = 10 ^ 

B. Density at Z 
From Table B- 

10810P ■ 
or 

p = 10“^ 


1.056 , 3 

gm/cm 


= 400 km 

5, for TE « 689 ,8 ^K, we have 
-15.285 

5.285/ , 3 

gm/cm 


At this height the seasonal-latitudinal correction for density 
according to Table B-4, is negligible. 
















RATIO OF THE LOCAL TEMPERATURE (TL) TO THE GLOBAL MINIMUM TEMPERATURE (TC) AS A 
FUNCTION OF LOCAL TIME AND OF LATITUDE (LAT). ALL RATIOS HAVE BEEN MULTIPLIED 
BY 1000 TO ELIMINATE THE DECIMAL POINT (AFTER REF. 1). 
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TABLE B-1. (continued) 
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TABLE B-1. (concluded) 
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TABLE B-2. 


TEMPERATURE INCREMENT AS A FUNCTION OF GEOMAGNETIC INDICES 

(AFTER REF. 1) 
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TABLE B-4. 


TABLES FOR THE SEASONAL-LATITUDINAL DENSITY VARIATION 
(AFTER REF. 1) ALOG p = (S)(P) sin2 (LAT) 


TABLE OF THE MAXIMUM AMPLITUDE, S » 0.02 (Z - 90) exp [-0.045 (Z - 901] 


Z(km) 

S 

Z (km) 

S 

Z (km) 

s 

90 

0. 000 

130 

B9 


0. 016 

95 

.0. 080 

140 



0. 007 

100 

0. 128 

150 

0. 081 

240 

0. 004 

105 

0. 153 

160 

0. 060 

260 


110 

0. 163 

170 

0. 044 

280 


115 


180 

0. 031 

300 

0. 000 

120 


190 

0. 022 




TABLE OF THE PHASE. P = sin 3^ (d + 100)* 


DATE 

P 

DATE 

P 

DATE 

P 

DATE 

P 

Jan. 1 

±0. 989 

Apr. 1 

TO. 129 

June 30 

TO. 994 

Sept. 28 

±0. 086 

1 1 

±0. 948 

1 1 

TO. 297 

July 10 

TO. 961 

Oct. 8 

±0. 255 

21 

±0. 880 

21 

TO. 456 

20 

TO- 900 

18 

±0. 417 

31 

±0. 786 

May 1 

TO. 602 

30 

TO. 812 

20 

±0. 567 

Feb. 10 

±0. 668 

11 

TO. 730 

Aug. ‘9 

TO. 699 

Nov. 7 

±0. 699 

20 

±0. 531 

21 

TO. 836 

19 

TO. 567 

17 

±0. 812 

Mar. 2 

±0. 378 

31 

TO. 918 

29 

TO. 417 

27 

±0. 900 

12 

±0. 214 

June 10 

TO. 972 

Sept. 8 

TO. 255 

Dec. 7 

±0. 961 

22 

±0. 043 

20 

TO. 998 

18 

TO. 086 

17 

±0. 994 

Apr. 1 

TO. 129 

30 

TO. 994 

28 

±0. 086 

27 

±0. 998 


•Take the upper sign for the Northern Hemisphere, the lower for the Southern Hemisphere. 


TABLE OF sin^ [LAT) 


LAT 

sin^ (LAT) 

LAT 

sin^ (LAT) 

LAT 

sin^ (LAT) 

0* 

0. 000 

o 

0 

0.250 

60® 

0. 750 

5 

0. 008 

35 

0. 329 

65 

0. 821 

10 

0- 030 

40 

0 413 

70 

0. 883 

15 

0. 067 

45 

0. 500 

75 

0. 933 

20 

0. 117 

50 

0. 587 

80 

0. 970 

25 

0. 179 

55 

0.671 

85 

0. 992 

30 

0. 250 

60 

0. 750 

90 

1. 000 
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TABLE B-5. 

ATMOSPHERIC DENSITY AS A FUNCTION OF HEIGHT AND EXOSPHERIC TEMPERATURE (DECIMAL 

LOGARITHMS, g/cm^) (AFTER REF. 1) 
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TABLE B-5. (continued) 
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APPENDIX C 
GLOSSARY* 


Atmospheric Density - Same as Mass Density. 

Celestial Longitude — The arc of the ecliptic between the vernal equinox and the point at which the celes- 
tial longitude is given. It is always measured eastward from the vernal equinox, completely around the 
ecliptic, from to 360°. 

Daily 10.7 cm Flux - Assumed to be the same as radio flux. A measured indicator for the amount of EUV 
solar radiation received by the Earth. See Radio Flux. 

Declination — In the geocentric coordinate system, it is the angular distance along the meridian of a point 
of a body from the equator. Declination is analogous to latitude on Earth. It is taken positive north of 
the equator and negative south of the equator. 

Density - See Mass Density and Number Density. 

Density Bulge — A slight bulge in the daylight portion of the atmosphere that is caused by atmospheric 
heating. The center of the bulge follows the Sun, lagging by two hours, and also migrates north and 
south with the sub-solar point. A slight depression in the dark portion of the atmosphere (anti-bulge), 
which is a product of the bulge, is centered eleven hours earlier at 0300 local time. At any given height 
above 120 km, the maximum density occurs at the center of the density bulge. 

Density Scale Height — The scale height of a point in the atmosphere is a numerical quantity that represents 
the altitude above the point at which the mass density would decrease by a factor of 1/e (the exponen- 
tial log e) from the density at the point. See equation A=35. 

Diffusive Equilibrium — The steady state resulting from the diffusive process in which the constituent 
gases of the atmosphere are distributed independently of one another. In such a state, the number 
density of the heavier constituents decreases more rapidly with altitude than that of the lighter constituents. 

Diurnal Effect - The day-to-night variation in nearly all atmospheric parameters that is caused by the rota- 
tion of the Earth. See Density Bulge. 

Ecliptic ~ The apparent path of the Sun about the Earth during a year. Strictly, it is the projection of the 
plane of the Earth’s orbit on the celestial sphere. 

Electromagnetic Radiation - Energy that is propagated through space, primarily from the Sun, in the form 
of an advancing disturbance in electric and magnetic fields (often called radiation). 

Exosphere - The outermost, or topmost, portion of the atmosphere. Its lower boundary is the critical 
level of escape, variously estimated at 350 to 750 km above the Earth’s surface. 


♦Cross references within the glossary are indicated by bold face. 
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Exospheric Temperature - The value that the atmospheric temperature reaches near 350 km altitude. 
Above this altitude, temperature is considered to be isothermal (lapse rate of zero) and to range from 
650° to 2100°K 

Geomagnetic Index - A measurement of the most active component of the magnetic field made by surface 
magnetic observatories (an average of 12 selected stations) at three-hour intervals. The measurement is 
made in the units of gauss or gamma (1 gauss = 10^ gammas), but observations are reported in terms of 
a unitless quantity (ap) which varies from 0 to 400. 

Gravity Wave — A wave disturbance in which buoyancy acts as the restoring force on parcels displaced from 
hydrostatic equilibrium. 

Heat Conductivity ~ An intrinsic physical property of a substance, describing its ability to conduct heat as 
a consequence of molecular motion. 

Heterosphere — The upper portion of the two-part division of the atmosphere according to the general ho- 
mogeniety of atmospheric composition; the layer above the homosphere. The heterosphere is character- 
ized by variation in composition, and mean molecular weight of constituent gases. The region starts at 
80 to 100 km above the Earth, and therefore closely coincides with the ionosphere and the thermosphere. 

Homosphere — The lower portion of a two-part division of the atmosphere according to the general homo- 
geniety of atmospheric composition; opposed to the heterosphere. The homosphere is the region in 
which there is no gross change in atmospheric composition, that is, all of the atmosphere from the 
Earth’s surface to about 80 to 100 km. 

Homopause - The top of the homosphere, or the level of transition between it and the heterosphere. It 
probably lies between 80 and 90 km, where molecular oxygen begins to dissociate into atomic oxygen. 
The homopause is somewhat lower in the daytime than at night. 

Hour Angle - The angular distance measured eastward or westward along the celestial equator to the longi- 
tude of the point for which the hour angle refers. Morning hour angles are negative and afternoon are 
positive. 

Hydrostatic Equilibrium — The state of a fluid whose surfaces of constant pressure and constant mass (or 
density) coincide and are horizontal throughout; a complete balance exists between the force of gravity 
and the pressure force. 

Ionization - In atmospheric electricity, the process by which neutral atmospheric molecules or other sus- 
pended particles are rendered electrically charged chiefly by collisions with high-energy particles. Cos- 
mic rays and emanations from radioactive gasses are the main sources of atmospheric ionization. In the 
lower atmosphere, decay electrons of mu-mesons plus alpha particles from radioactive gases, as well as 
beta particles and gamma rays, serve to ionize air molecules. 

Julian Date (also Julian day) — Number of days measured from January 1 (noon GMT), 4713 B.C. 

Lapse Rate - The change in temperature with increasing altitude. 

Magnetic Field — A region wherein any magnetic dipole would experience a magnetic force or torque. 
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j Mass Density — The ratio of the mass of any substance to the volume occupied by it (usually expressed in 

1 gm/cm^). 

I 


Mean Free Path — The average distance traveled by the molecules of a perfect gas between consecutive colli- 
sions with one another. 

Mesosphere — The atmospheric shell between about 20 km and about 70 or 80 km, extending from the top 
of the stratosphere to the upper temperature minimum (the mesopause). It is characterized by a broad 
temperature maximum (the mesopeak) at about 50 km, except possibly over the winter polar regions. 

Mixing — A random exchange of atmospheric constituents caused primarily by non-homogene ous pressure 
forces. 

Modified Julian Day - Number of days measured from November 17, 1858 (midnight GMT). 


Number Density — The numerical count of molecules of a particular constituent for a given volume (usually 
given as number/cm^). 


Photodissociation — The dissociation (splitting) of a molecule by the absorption of a photon. The resulting 
components may be ionized in the process (photoionization). 

Radio Flux (F10.7) — The radio flux density at 10.7 cm is a useful indicator of solar activity as it exhibits 
both 11-year and 27-day periodicities. Its high correlation with the solar radiation which is absorbed 
in the upper atmosphere makes it a desirable weighting function for representing the effects of this ab- 
sorption in a model atmosphere. The standard data source is Ottawa, Canada, although it is measured at 
several observatories. The solar flux is actually measured over a complete bandwidth to increase the 
faithfulness of the radio energy input as it passes through the receiver and therefore must be divided by 
width of the band. The 10.7 cm solar flux has units of lO"^^ watts/sq m/ sec/bandwidth but should be 
considered as a unitless quantity for the equations herein. 

Right Ascension — The angular distance measured from the vernal equinox eastward along the celestial 
equator to the longitude of the point to which the right ascension refers. 

Solar Radiation — The total electromagnetic radiation and corpuscular radiation emitted by the Sun. 


Solar Wind - The flux of plasma from the Sun. 
Static Diffusion — Same as Diffusive Equilibrium. 


Stratosphere — The atmospheric shell above the troposphere and below the mesosphere. It extends, there- 
fore, from the tropopause to the height where the temperature begins to increase in the 20 to 25 km 
region. 


i 
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Sunspot Number — A solar index which has been compiled back to 1600 A.D. The sunspot number takes 
into account the number of sunspot groups as well as the number of individual spots. It can be com- 
puted from the formula 


R = K(10g + s) 


where 


R = sunspot number 
g = number of groups 
s = number of spots 
K = a constant, roughly equal to 1 

K is used to adjust individual R numbers to account for some observatories having better observing con- 
ditions than others. 

Thermosphere - The region of the atmosphere extending upward from about 85 km to the exosphere. In 
this region the temperature increases with altitude to about 350 km and then is constant with increasing 
altitude. 

Thermopause — The top of the thermosphere, or the level of transition between it and the exosphere. 

81-Day Mean Solar Flux - The arithmetic average of the daily 10.7 cm solar flux values for the 81 days 
preceding the day for which the 81-day mean is given. 
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NASA SPACE VEHICLE DESIGN CRITERIA 
MONOGRAPHS 


ENVIRONMENT 

SP-8005 

SP-8010 

SP-8011 

SP-8013 

SP-8017 

SP-8020 

SP-8021 

SP-8023 

SP-8037 

SP-8038 

SP-8049 

SP-8067 

SP-8069 

SP-8084 

SP-8085 

SP-8091 

SP-8092 


Solar Electromagnetic Radiation, revised May 1971 

Models of Mars Atmosphere (1967), May 1968 

Models of Venus Atmosphere (1972), revised September 1972 

Meteoroid Environment Model- 1969 (Near Earth to Lunar Sur- 
face), March 1969 

Magnetic Fields— Earth and Extraterrestrial, March 1969 
Mars Surface Models (1968), May 1969 

Models of Earth’s Atmosphere (190 to 2500 km), revised March 1973 
Lunar Surface Models, May 1969 

Assessment and Control of Spacecraft Magnetic Fields, September 
1970 

Meteoroid Environment Model— 1970 (Interplanetary and Plan- 
etary), October 1970 

The Earth’s Ionosphere, March 1971 

Earth Albedo and Emitted Radiation, July 1971 

The Planet Jupiter (1970), December 1971 

Surface Atmospheric Extremes (Launch and Transportation Areas), 
May 1972 

The Planet Mercury (1971), March 1972 
The Planet Saturn (1970), June 1972 

Assessment and Control of Spacecraft Electromagnetic Inter- 
ference, June 1972 
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SP-8103 

SP-8105 

STRUCTURES 

SP-8001 

SP-8002 

SP-8003 

SP-8004 

SP-8006 

SP-8007 

SP-8008 

SP-8009 

SP-8012 

SP-8014 

SP-8019 

SP-8022 

SP-8029 

SP-8031 

SP-8032 

SP-8035 

SP-8040 

SP-8042 

SP-8043 


The Planets Uranus, Nepture, and Pluto (1971), November 1972 
Spacecraft Thermal Control, May 1973 

Buffeting During Atmospheric Ascent, revised November 1970 

Flight-Loads Measurements During Launch and Exit, revised June 
1972 

Flutter, Buzz, and Divergence, July 1964 
Panel Flutter, revised June 1972 

Local Steady Aerodynamic Loads During Launch and Exit, May 
1965 

Buckling of Thin-Walled Circular Cylinders, revised August 1968 

Prelaunch Ground Wind Loads, November 1965 

Propellant Slosh Loads, August 1968 

Natural Vibration Modal Analysis, September 1968 

Entry Thermal Protection, August 1968 

Buckling of Thin- Walled Truncated Cones, September 1968 

Staging Loads, February 1969 

Aerodynamic and Rocket-Exhaust Heating During Launch and 
Ascent, May 1969 

Slosh Suppression, May 1969 

Buckling of Thin-Walled Doubly Curved Shells, August 1969 
Wind Loads During Ascent, June 1970 
Fracture Control of Metallic Pressure Vessels, May 1970 
Meteoroid Damage Assessment, May 1970 
Design— Development testing. May 1970 
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SP-8044 

SP-8045 

SP-8046 

SP-8050 

SP-8053 

SP-8054 

SP-8055 

SP-8056 

SP-8057 

SP-8060 

SP-8061 

SP-8062 

SP-8063 

SP-8066 

SP-8068 

SP-8072 

SP-8077 

SP-8079 

SP-8082 

SP-8083 

SP-8095 


Qualification testing, May 1970 
Acceptance testing, April 1970 

Landing Impact Attenuation for Non-Surface-Planing Landers, 
April 1970 

Structural Vibration Prediction, June 1970 

Nuclear and Space Radiation Effects on Materials, June 1970 

Space Radiation Protection, June 1970 

Prevention of Coupled Structure-Propulsion Instability (Pogo), 
October 1970 

Flight Separation Mechanisms, October 1970 

Structural Design Criteria Applicable to a Space Shuttle, revised 
March 1972 

Compartment Venting, November 1970 

Interaction with Umblicals and Launch Stand, August 1970 

Entry Gasdynamic Heating, January 1971 

Lubrication, Friction, and Wear, June 1971 

Deployable Aerodynamic Deceleration Systems, June 1971 

Buckling Strength of Structural Plates, June 1971 

Acoustic Loads Generated by the Propulsion System, June 1971 

Transportation and Handling Loads, September 1971 

Structural Interaction with Control Systems, November 1971 

Stress-Corrosion Cracking in Metals, August 1971 

Discontinuity in Metallic Pressure Vessels, November 1971 

Preliminary Criteria for the Fracture Control of Space Shuttle 
Structures, June 1971 
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SP-8099 Combining Ascent Loads, May 1972 

GUIDANCE AND CONTROL 

SP-8015 Guidance and Navigation for Entry Vehicles, November 1968 

SP-8016 Effects of Structural Flexibility on Spacecraft Control Systems, 

April 1969 

SP-8018 Spacecraft Magnetic Torques, March 1969 

SP-8024 Spacecraft Gravitational Torques, May 1969 

SP-8026 Spacecraft Star Trackers, July 1970 

SP-8027 Spacecraft Radiation Torques, October 1969 

SP-8028 Entry Vehicle Control, November 1969 

SP-8033 Spacecraft Earth Horizon Sensors, December 1969 

SP-8034 Spacecraft Mass Expulsion Torques, December 1969 

SP-8036 Effects of Structural Flexibility on Launch Vehicle Control Sys- 

tems, February 1970 

SP-8047 Spacecraft Sun Sensors, June 1970 

SP-8058 Spacecraft Aerodynamic Torques, January 1971 

SP-8059 Spacecraft Attitude Control During Thrusting Maneuvers, February 

1971 

SP-8065 Tubular Spacecraft Booms (Extendible, Reel Stored), February 

1971 

SP-8070 Spaceborne Digital Computer Systems, March 1971 

SP-8071 Passive Gravity-Gradient Libration Dampers, February 1971 

SP-8074 Spacecraft Solar Cell Arrays, May 1971 

SP-8078 Spaceborne Electronic Imaging Systems, June 1971 

SP-8086 Space Vehicle Displays Design Criteria, March 1972 
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SP-8096 


Space Vehicle Gyroscope Sensor Applications, October 1972 


SP-8098 Effects of Structural Flexibility on Entry Vehicle Control Systems, 

June 1972 

SP-8102 Space Vehicle Accelerometer Applications, December 1972 

CHEMICAL PROPULSION 

SP-8025 Solid Rocket Motor Metal Cases, April 1970 

SP-8039 Solid Rocket Motor Performance Analysis and Prediction, May 

1971 

SP-8041 Captive-Fired Testing of Solid Rocket Motors, March 1971 

SP-8048 Liquid Rocket Engine Turbopump Bearings, March 1971 

SP-805 1 Solid Rocket Motor Igniters, March 1 97 1 

SP-8052 Liquid Rocket Engine Turbopump Inducers, May 1971 

SP-8064 Solid Propellant Selection and Characterization, June 1971 

SP-8075 Solid Propellant Processing Factors in Rocket Motor Design, 

October 1971 

SP-8076 Solid Propellant Grain Design and Internal Ballistics, March 1972 
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NASA-Langley, 1973 


